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Abstract. Numerical aspects of dynamos in periodic domains are discussed. Modifications of 
the solutions by numerically motivated alterations of the equations are being reviewed using 
'f0^', the examples of magnetic hyperdiffusion and artificial diffusion when advancing the magnetic 
I . field in its Euler potential representation. The importance of using integral kernel formulations 
^~>) ' in mean-field dynamo theory is emphasized in cases where the dynamo growth rate becomes 
, comparable with the inverse turnover time. Finally, the significance of microscopic magnetic 
^ Prandtl number in controlling the conversion from kinetic to magnetic energy is highlighted. 
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^ 1. Introduction 

I— I There are two important aspects connected with astrophysical dynamos compared with 

\ dynamos on a bicycle. Firstly, they are self-excited and do not require any permanent 
^0 ■ magnets. Secondly, they are homogeneous in the sense that the medium is conducting 
' everywhere in the dynamo proper and there are no wires or insulators inside. Self-excited 
CIh dynamos were invented by the Danish inventor S0ren Hjorth, who received the patent 
Q for this discovery in 1854, some 12 years before Samuel Alfred Varley, Ernst Werner 
$^ von Siemens and Charles Wheatstone announced such an invention independently of 
c/3 each other. Von Siemens is known for having recognized its industrial importance in 
^ ^ ^ producing the most powerful generators at the time, for which he, in turn, received a 
patent in 1877. 

The idea th at homogeneous dynamos might work in the Sun , was first prop osed by 
Larmoij (|l919l ) in a one-page paper. However, some 14 years later. ICowlinj (119331 ) showed 
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that axisymmetric dynamos cannot work in a body like the Sun. At the time it was 
not clear whether this failure was genuine, or whether it was critically connected with 
Cowling's assumption of axisymmetry. The sus picion th at the third dimension might be 



critical was not particularly emphasized when iLarmoij (|1934I ) tried to defend his early 
I . suggestion with the words "the self-exciting dynamo analogy is still, so far as I know, 
I the only foundation on which a gaseous body such as the Sun could possess a magnetic 
I . field: so that if it is demolished there could be no explanation of the Sun's magnetic field 
' even remotely in sight. 



The essential idea about the operation of the solar dynamo came from |Parkeil(|l955l) . 



' who developed the notion that cyclonic events would tilt a toroidal field systematically 
^ . in the poloidal direction, closing thereby a critical step in the d ynamo cycle . While this 
■ " " ' concept is still valid todav. it still required the existence proof bv lHerzenber j ( 19581 ) that 



began to convince critics that Cowling's antidynamo theorem does not extend to the 
general case of three dimensions. 

Nevertheless, subsequent progress in modeling the solar dynamo appears to have been 
suspended un til the foundations of a mean-field treatment of the induction equation were 
developed by ISteenbeck et al.l ( 1966 ) . In the following year s, a large number of models 



were computed covering mostly aspects of the solar dynamo ( Steenbeck fc Krauselll969a : 
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Parkeill970allbll3.ll971blldl l?). but in some cases also terrestrial d ynamos (ISteenbeck &: Krausd 
1969bHParkeilll97ld) and the galactic dynamo (|Par ke3ll971 aUel: IVamshtein &: RuzmaikinI 



1971 LTi972|) . These developments provided a major boost to dynamo theory given that 



until then work on the gal actic dynamo , for example, focussed on aspects concerning the 
small-scale magnetic field ( Parker 19691 ) , but not the global large-scale fields on the scale 
of the entire galaxy. In fact, also regardi ng small-scale dynamos, there were important 
developments made by iKazantsev ( 1968 ) , but t hey remained mo s tly un noticed in the 
West, even when the first direct simulations bv iMeneguzzi et"aL (1981) demonstrated 
the operation of such a dynamo in some detail. In fact, in some of these dynamos, the 
driving of the flow involved helicity, but its role in helping the dynamo remained un- 
convincing, because no large-scale field was produced. We now understand that this was 
mainly because there was not enough scale separation betwee n the scale of the d omain 
and the forcing scale, and that one needs at least a ratio of 3 (jHaugen et al.ll2004l ). 

Simulations in spherical geometry we re much more readily able to demonstrate the 
production of large-scale magnetic fields ( Gilman 1983 : Glatzmaieilll985 ). b ut even today 



these simulations produce magnetic fields that propagate toward the poles (jKapyla et al 



20ld) and not toward the equator, as in the Sun. We can only speculate about possible 



shortcomings of efforts such as these that must ultimately be able to reproduce the solar 
cycle. 

Several important developments happened in the 1980s. Firstly, it be came broadl y 
accepted that the magnetic field inside the Sun might be in a fibril state (Parker|ll98i), 
i.e. the filling factor is small and most of the field is concentrated into thin flux tubes, 
as manifested by the magnetic fleld appearance in the form of sunspots at the surface. 
However, such tubes would be magneti cally buoyant, a.nd are expe cted to rise to the 
surface on a time scale of some 50 days (|Moreno-Insertislll983l Il986f) . This time is short 
compared with the cycle time and might lead to excessive magnetic flux losses, which then 
led to the proposal that the magnetic field would instead be generated in the overshoot 
layer beneath the convection zone. This idea is still the basic picture today, although 
simulations of convection generally produce magnetic fields that are distributed over the 
entire convection zone. 

Yet another important development in the 1980s was the proposal that the a effect 
might actually be the sum of a kinetic and a magnetic part and that the magnetic part 
can be estimated by solving an evolution equation for the magnetic helicity density. The 
importance of this development was obscured by the excitement that the two evolution 
equations for poloidal and toroidal field, sup plemented by a t hird equation for the mag- 
netic helicity density, could produce chaos (jRuzmaikinI 1 1 98 II ) . The connection to what 
was to come some 10-20 years later was not yet understood at that point. Simulations 
of helical MHD turbulence in a periodic domain demonstrated that in a periodic domain 
the a effect might be quenched in an Rej\/-dependent fashion like 
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(1.1) 



If this were also true of astrophysical dynamos, a would be negligibly small and would not 
be relevant for explaining the magnetic field in these bodies. Such quenching is therefore 
nowadays referred to as catastrophic quenching. However, there is now mounting evidence 
that this type of a quenching is a special case of a more general formula (see, e.g., 
lBrandenbur3l2008h 
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RcAf [mMo J • BjBl^ - (V • Fc)/2kfBl^ ~ {da/dt)/2r]tkf] 



RCM-B2/B2 



(1.2) 
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which comes from magnetic hehcity conservation. Note that in this equation there are 3 
new terms that ah scale with ReM and are therefore important. Even in a closed or peri- 
odic domain, the first and third terms in squared brackets contribute, the most promising 
way out of catastrophic qu enching is through magnetic helicity fluxes ( Blackman fc Fieldl 
2OOOI: iKleeorin et al. I l2000l). These developments are still ongoing and we refer here to 



some r ecent papers by Shukurov et al. ( 2006 ): Brandenburg et al. ( 20091 ): Candelaresi et al, 
( 2010h . 



2. Simulating dynamos 

2.1. Roberts flow 

Much of the theoretical understanding of dynamos is being helped by numerical simu- 
lations. In fact, nowadays one of the simplest dynamos to simulate is the Roberts flow 
dynamo. A kinematic dynamo can be simulated by adopting a velocity field of the form 

U = V X ijjz + ki'4)z, (2.1) 

where ip = {Uq / ko) cos kox cos koy . In order to give some idea about the ease at which 
reasonably accurate solutions can be obtained we give in Table[l]the numerically obtained 
critical values of the magnetic Reynolds number. Re a/ = itj-ms/'y^f, for a given resolution. 

Table 1. Critical values of Rbm for the Roberts flow dynamo at low resolution (from 8'^ to 
32^ mesh points) and difi'erent spatial order of the numerical scheme. 



Resolution 



Re 



A/,crit 



2nd order 6th order 8th order 



8^ 
32^ 



5.23 5.15 5.16 

5.517 5.514 5.518 
5.522 5.522 5.521 



Even turbulent dynamos are nowadays easy to simulate and meaningful results have 



been obtained already at relatively low resolution, provided the flow is helical (jBrandenbure 



20011 ). However, there are also examples where numerical aspects can have a major effect 
on the outcome of such simulations. In the following we discuss two examples: magnetic 
hyperdiffusion and the use of Euler potential with artificial magnetic diffusion. 

2.2. Magnetic hyperdiffusion in helicity-driven dynamos 

Dynamos work by maintaining the magnetic field against Ohmic decay via magnetic 
induction. It is then not surprising that the result can be sensitive to the numerical 
treatment of magnetic diffusion. One example is the consideration of magnetic hyperdif- 
fusion of the form 

77J ^ (-l)"-ii/„V2"-2 J (2.2) 

instead of the regular rjj term. In helical dynamos in closed or periodic domains the 
saturation time is given by rjk^At = 1, but with hyperdiffusion this condition becomes 
r]nki"At — 1; see Figure [T] This is exactly what one expects from a scheme like hyperdif- 
fusion that enhances the effective magnetic Reynolds number. H owever, another impor 



tant m odification is that the saturation amplitude increases from ([Brandenburg fc Sarson 



2002h 



^.ftV". (2.3) 
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Figure 1. Evolution of large-scale magnetic field (dashed lines) for runs with regular mag- 
netic diffusion (n = 1, left, «rnis/f?fci ~ 45) and magnetic hyperdiffusion (n = 3, right, 
Ut ths I Vzki = 3200), at 16^ resolu tion and fcf/fci = 3, compared with the prediction (solid lines) 
of [Brandenburg fc SarsonI (|2002f ). where ts is the saturation time of the small-scale field. 



This is a caveat that is important to keep in mind when employing magnetic hyperdiffu- 
sion for astrophysical simulations. This has a major effect on the saturation of large-scale 
magnetic fields that is at first glance surprising. However, once one realizes that the sat- 
uration in a periodic domain is governed by the magnetic helicity equation applied to 
the steady state, i.e. by 

^(A.B) = -277(J.B), (2.4) 

where we split the right-hand side into contributions from large- and small-scale fields, 
i.e. B = B + b and J = J + j, so that {J ■ B) — {J ■ B) + {j -b), we have for fully helical 
large- and small-scale fields in the steady state, 

= -277(7 • - 27j{j ■ b) « ±277 (fci (B') - fcf (b')) , (2.5) 

with (b^) « i?^q, it becomes clear that the use of magnetic hyperdiffusion picks up the 
fci and ki factors at correspondingly higher powers, leading thus to Equation (|2.3p . The 
upper and lower signs of the term on the right-hand side of Equation (12. 5p apply to 
small-scale forcings with positive and negative helicity, respectively. 

2.3. MHD with Euler potentials 
Until recently, the use of Euler potentials (EP) has been a pop ular choice for solv- 



ing the MHD equatio ns numerically using Lagrangian methods ([Price fc Bate! 12007 



Rosswog fc Pricdl2007l) . The representation of B in terms of EP, a and /?, as 

B = VaxVI3 (2.6) 

is a nonlinear one, which solves the induction equation in the case 77 = 0, dB/dt = 
^ X {U X B), provided Da/Dt = D/?/Dt — 0. The problem is that the magnetic field 
tends to develop sharp structures that will not be properly resolved by the numerical 
scheme. The hope has been that the overall properties of the magnetic field at larger 
scales would then still be approximately correct. However, this turns out not to be the 
case. And, more importantly, as one increases the resolution, the solution converges, but 
it is simply the wro ng solution. This has been demonstrated in detail in a separate paper 



( Brandenburg! [2OIO ) where, among other cases, solutions of the Galloway-Proctor fast 



dynamo flow were considered. In Figure [2] we show that the solution of the induction 
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Figure 2. Left: Comparison of the evolution of Brms for the modified GaUoway-Proctor flow 
with point-wise zero hehcity for methods A and EP using 256"^ meshpoints and ReM = 10*. 
Note the power law scaling for the EP method and the exponential scaling for the A method. 
Right: ReAf dependence of the exponents A and a characterizing the evolution of Brms ~ e^' 
for the A method and Brms ~ t" for the EP method for the modified Galloway-Proctor flow 
with point-w i se zer o helicity for methods A and EP using 256"^ meshpoints. Adapted from 
lBrandenbur"3 l|2010l ). 



equation using EP leads to algebraically decaying solutions of the form 



Br 



t" (EP method), 



(2.7) 



where, with increasing resolution, a converges toward a value around —3, while with 
the usual vector potential method (A method) one solves dA/dt — U x B — rjj with 
B — V X A and J — ^ x B, and finds instead 



Bin 



(A method). 



(2.8) 



where A also converges, and its value is about 0.22 for Rcm — lO'*; see Figure[2] 

It is quite remarkable that by changing the properties of the solution at small scales 
only slightly, one can produce rather dramatic effects. This includes cases of magnetic 
hyperdiffusivity, where the large-scale field amplitude can be quite different, albeit in 
agree ment with the theory applied to the hyperdifTusive case ([Brandenburg fc Sarson 
[2003). Another example is that of artificial diffusion in solutions for the EP, where the 
obtained results bear no resemblance with those obtained using the A method. 

2.4. Quantitative comparison between simulations and mean-field theory 

Mean-field theory has the potential of being a quantitatively accurate and hence predic- 
tive theory. In order to establish this in particular cases, it is important to consider as 
many contact points between theory and simulations as possible. One thing we can do is 
to determine a effect, turbulent diffusion, and other effects from simulations, and to com- 
pare the thus calibrated mean-field model with simulations. This provides an important 
resource for ideas of what one might have been missing in various contexts. Here we just 
mention the case of the Roberts flow, for whi ch aij and rjij have been determined using 



the test-field methods (jSchrinner et al.l 120071 ) . One might then expect that the growth 



rate obtained from the underlying dissipation rate, 

A = afc — (77t + ?y)fc^ 



(2.9) 



should agree with the value obtained from the direct calculation. This is however only 
the case for A = 0, but not for A 7^ 0. However, it would be a mistake to assume that 
there is something wrong with the test-field methods. Instead, what is wrong here is just 
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Figure 3. Left: Hsm dependence of the growth rate for the Roberts flow as obtained 
from a direct calculation (Agrowth) compared with the result of the dispersion relation, 
Adisp = akz — {ri + r]t)k^, using a cubic domain of size L'^, where fci = 2-k/L and fcf = \/^ki. For 
this range of Rcm, the most unstable mode is the largest one that fits in the box {kz — ki). The 



two horizontal lines in gray mark the values of Ag 



and Adisp at ReAf = 30, denoted by (i) and 



(ii), respectively. Right: Laplace-transformed effective growth rate, A(s) — a{s)k— [rj + rit{s)]k'^ , 
for the Roberts flow with ReM = 10 and 50. Note the different signs of the slope at the inter- 
section with the diagonal (denoted by circles). Adapted from Hubbard & Brandenburg (200^). 

the assumption of homogeneity and stationarity. Obviously, when A ^ 0, the solution is 
exponentially growing or decaying like e^*, so it is clearly not steady! 

When the assumption of steadiness is no longer satisfied, one has to return to the 
underlying integral relation between a and the mean field, i.e. a{z,t)B{z,t) has to be 
replaced by a convolution, so 

a{z,t)B{z,t) ^ J a{z- z',t-t')B{z',t')dz'dt', (2.10) 

and similarly for the magnetic diffusivity. One way of dealing with this complication is to 
note that the convolution in real space corresponds to a multiplication in Fourier space. 
In other words, we can write 

£{k,io) ^ a{k,uj)B{k,io), (2.11) 

where £{k,u>) = J 5(z, i) e"**^'^^^'^*Mz dt is the Fourier-transformation of £{z,t) (and 
likewise for the other fields). Of course, the value of w that is of interest is = iA, where 
A is the then self-consistently obta ined growth rate. This has been described in detail 
bv lHubbard fc Brandenburg (|20Q9t ). who motivated their study using the example of the 
Roberts flow, where the discrepancy between the numerical solution and that obtained 
for A = is shown. A reasonable fit to their data is a{k,Lo) = 1/(1 — iwr) for k ^ 0. 



3. Low magnetic Prandtl number and application to accretion discs 

In many astrophysical bodies the magnetic Prandtl number, Pi'm ~ v/fj, is either large 
or small, but not around unity. Again, from a numerical point of view, it is surprising 
that the ratio v/rj is important even though 

v^O, rj^O. (3.1) 

In many numerical treatments it is implicitly assumed that the exact values of v and ry 
do not explicitly matter, because it should not matter how long each of the two turbulent 
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cascades is. This should indeed be true provided the dynamics we are interested in takes 
place entirely on the large scales. But for helical large-scale dynamos, this is evidently 
not the case, because in the kinematic regime, the magnetic energy spectrum follows 
the Kazantsev k^^"^ spectrum and reaches a peak at the resistive scale, corresponding to 
the wavenumber = (Jmis/^)- However, this applies first of all only to the kinematic 
regime and, seemingly, it is relevant only for the small-scale dynamo regime. This has been 
discussed recently in connection with the contrasting case of l arge-scale dynamos , wher e 
the excitation conditions are not affected by the value of Prjvf ( Brandenburgll2009l 2011 ). 
Nevertheless, even in this case the dissipation rates of kinetic and magnetic energies, 
CK and eM, respectively, do depend on the value of Ptm- Although it seems fairly clear 
that the ratio ex/fM increases with Ptm in power law fashion proportional to PrJ^, the 
exponent n is not well constrained. Earlier work covering the range 10~^ ^ PrM ^ 1 
at 512"^ resolution suggested n — 1/2, although additional data covering also the range 
1 < Pi'M ^ 10'^ given an exponent closer to n = 0.6 or even n = 2/3. There is at present 
no theory for the value of the exponent n. 

Applying these finding s to accretion discs, we should first recall that the work of 



Lesur fc Longarettil (|2007[ ) suggested that the onset of the magneto-rotational instability 



shows a strong dependence on PrM- However, one should expect that when large-scale 
dynamo action is possible, this conditi on may change and the onset would then be in- 
dependent of PrM ■ This is indeed what iKapvla fc Koroil tOKt find. The question is of 
course what is the relevant large-scale dynamo mechanism in this case. One proposal is 
the incoherent a-shear dynamo, which works through constructive amplification of the 
mean field in the direction of mean shear. Yet an other possible mechanism is the shear- 
current dynamo ( Rogachevskii fc Kleeorin 2003f) . although this result has not yet been 
confirmed. 



4. Conclusions 



Simulating simple dynamos on the computer is nowadays quite simple. Nevertheless, we 
have seen here examples that illustrate that things can also go quite "wrong". In the case 
of magnetic hyperdiffusion, it is clear what happens (jBrandenburg fc SarsonI 120021 ) . so 
th at magnetic hyperdi f fusion can also be used to ones advantage, as was demonstrated 
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()2002l) . However, in the case of Euler potentials it is not clear 
what happens and whether this method can be used to simulate even the ideal MHD 
equations, given that each numerical scheme will introduce some type of diffusion. In this 
short review, we have also attempted to clarify why nume rical calculations of a effect 
and turbulent diffusion using the standard test-field method ( Schrinner et al ]|2007l) would 
yield values that can only reproduce a correct growth rate in the case of vanishing growth. 
In all other case, a representation in terms of integral kernels has to be used. Finally, we 
have discussed some effects of using magnetic Prandtl numbers that are different from 
unity. It turns out that in the steady state, the rate of transfer from kinetic to magnetic 
energy depends on the value of Pr^f- This is somewha t unexpected, beca use the onset 
condition for dynamo action docs not depend on Pr^f ( Brandenbur j 2009f) . and yet the 
actual efficiency of the dynamo, as characterized by the work done against the Lorentz 
force, —{U ■ (J X B)), does depend on Pr^f and is proportional to Pr^" (with n between 
1/2 and 2/3) for large values of Pr 

Understanding the limits of numerical simulations is just as important as appreciating 
its powers. As the example with the problem with mean-field and simulated growth 
rates shows, understanding the initial mismatch can be the key to a more advanced and 
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more accurate theory that will ultimately be needed when describing some of the yet 
unexplained properties of astrophysical dynamos. 
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